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^ . Abstract 

^ • Consider in ^^(R'^), rf > 1, the operator family H{g) := Hq + igW. Hq = ajai + 

. . . + ajfld + d/2 is the quantum harmonic oscillator with rational frequencies , W aV 
symmetric bounded potential, and g a real coupling constant. We show that if |f;| < p, 
p being an explicitly determined constant, the spectrum of H{g) is real and discrete. 
Moreover we show that the operator H{g) = a*ai + 0*20,2 + iga*20\ has real discrete 
spectrum but is not diagonalizable. 
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?H : 1 Introduction 

A basic fact underlying "PT-symmetric quantum mechanics (see e.g. [1-10]; V is the 

\ parity operation, and T the complex conjugation) is the existence of non self-adjoint, and 

not even normal, but "PT-symmetric Schrodinger operators (a particular case of complex 
CN ■ 

Tij" ' symmetric operators, as remarked in [11]) which have fully real spectrum. 

\ Two natural mathematical questions arising in this context are (i) the determination 

■ of conditions under which PT-symmetry actually yields real spectrum (for results in this 
direction see e.g. [12j.[13]. |14j . |15j . |16j ) and (ii) the examination of whether or not this 

^ , phenomenon can still be understood in terms of self-adjoint spectral theory; for example, 

■ it has been remarked that if a T^T-symmetric Schrodinger operator with real spectrum is 
diagonalizable, then it is conjugate to a self-adjoint operator through a similarity map (see 
e.g. [17] . |18j . [19j). Hence the question arises whether PT-symmetric Schrodinger- type 
operators with real spectrum are always diagonalizable. 

In this paper a contribution is given to both questions. First, we solve in the negative 
the second one. Namely, we give a very simple, explicit example of a VT symmetric 
operator, with purely real and discrete spectrum, which cannot be diagonalized because 
of occurrence of Jordan blocks. The example is the following Schrodinger operator, acting 
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in a domain D{P{g)) C L^(R^) to be specified later: 



H{g) ■.= a\ai + a*2a2 + iga*2ai + l, ge 



XI) 



Here ai, a* , i = 1,2 are the standard destruction and creation operators of two independent 
harmonic oscillators: 



1 

71 



Xi + 



d 

dxi 



1 

71 



so that can be rewritten under the form 



H{9) 



dxl 



+ xl 



dX2 



+ xl 



+ ^2 I ^2 



d 

dxi 



d 
dx2 



(1.2) 



Xl + 



:i.3) 



which is manifestly invariant under the "PT-operation X2 — > —X2, ig —ig- 

Second, we identify a new class of non self-adoint, "PT-symmetric operators with purely 
real spectrum in L^(R'^), d > 1. To our knowledge, this is the first example of such 
operators in dimension higher than one (a preliminary version of this result, without 
proofs, already appeared in [20]). An example of an operator belonging to this class is 
represented by a perturbation of the harmonic oscillators in dimension higher than one, 
namely by the following Schrodinger operators: 



1 



H{9) = ^Y. 



k=l I 



d^ 



dx 



I 2 2 
2 + ^k^k 



+ igW{xi,...,Xd) 



(1.4) 



Here W G L°°(E'^), W{-xi, . . . , -x^) = -W{xi, . . . ,Xd), \g\ < p, where p > is an 

explicitly estimated positive constant, and the frequencies uj^ > are rational multiples 

Pk 



of a fixed frequency cj > 0: i^k 



Qk 



"-L0. Here pk G N,qk G N : k 



1, . . . , d is a pair of 

UJl p 



relatively prime numbers, with both p^ and odd, k = 1, . . . ,n. When d = 2 

UJ2 q 

this result can be strenghtened: if uji/uj2 = p/q, the spectrum is real if and only if p and 
q are both odd. 

The paper is organized as follows: in the next section we work out the example (jl.ip 
making use of the Bargmann representation, in Section 3 we establish the class of "PT- 
symmetric operators with real spectrum by exploiting the real nature of Rayleigh-Schro- 
dinger perturbation theory (for related work on spectrum of PT-symmetric operators 
through perturbation theory, see [22], [23]), and in Section 4 we work out the example 
represented by the perturbation of the resonant harmonic oscillators proving the above 
statements. 



2 



2 A non diagonalizable VT symmetric operator with real 
discrete spectrum 

Consider the operator ii{g) whose action on its domain is specified by (II. ip or, equivalently, 
(jl.3p . Denote the operator corresponding to the two-dimensional harmonic oscillator, 
namely: 



1 


S 2' 


1 


2' 


2 




+ 2 





(2.11 



It is immediately verified that Vu := 0201^ G if u € D{Hq). Therefore we can give the 
following 

Definition 2.1 The operator family H{g) : g ^ M in L^(IR^) is the operator H{g) whose 
action is Hq + igV on the domain D{Hq). 

Then we have: 

Theorem 2.2 Consider the operator family H{g) defined above. Then, G M, j^j < 2: 

1. H{g) has discrete spectrum. 

2. All eigenvalues of H{g) are Am = m + l,m = 0, 1,2, . . .. Each eigenvalue Xm has 
geometric multiplicity 1 but algebraic multiplicity m + 1. 

More precisely: for each m there is an m- dimensional subspace Tim invariant under H{g) 

00 

such that we have the orthogonal decomposition = Tim! ^/ we denote Hm '■= ^^Iw™ 

m=0 

00 

the restriction of H{g) to Tim, then H{g) = Hm O'^d Hm is represented by the {m + 

m=0 

1) X (m +1) matrix: 



Hm = (m + l)/( 

m+l)x(m+l) 

+ igD 

n 

Here Dm is a nilpotent of order m + 1 . Explicitly: 



(2.2) 



/ ^ 

v/2(m - 1) 



V 



\ 




V3(m - 2) 



m 
/ 



D 



rn+l 



(2.3) 
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Remarks 

1. Spec{H{g)) is thus real and independent of g. 

2. Formula (j2.2p is the Jordan canonical form of Hm- The algebraic muliplicity is m+1. 
Since Dm 7^ 0, Hm is not diagonalizable by definition and, a fortiori, neither is H{g). 

Proof of Assertion 1 

The classical Hamiltonians corresponding to the operators Hq and H{g) represent their 
symbols, denoted ctqI^^jO ^-^d o"g(j;,^), respectively: 

^o(x,0 = + + A + (2-4) 

<^g{x,0 = ao{x,0 +'i9^{x,0, ^{x,0 ■= ^(^2 - iC2){xi + i^i) (2.5) 

We have indeed (formally) ao{x,—iVx) = Hq, ag{x,—iVx) = H{g). Since do +00 as 
1^1 + |x| +00, by well known results (see e.g.|24J, §XIII.14) it is enough to prove that 
^Isl < 5* = 2, and V (x, ^) outside some fixed ball centered in the origin of M^: 

0<il-^\g\)aoix,0 < Wg{x,0\ (2.6) 

To see this, we estimate: 

\^\ < ^\X2 -i(.2\\xi < ^i\x2 + \X1 +«Cl|^) = ^CTq, 

and hence 

kgl > |<To| - \g\\a\ > (1 - ■y)o'o- 

This proves the inequality and hence the assertion. 

To prove the remaining assertions of the theorem we make use of the Bargmann rep- 
resentation [21]. To this end, recall the general definition of the Bargmann transform Ub 
(even though we shall need it only for d = 2): 

{Ubu){z) := fiz) = j:^ f e-^'+2^<^.'?)-'?\(g) dq, z € & (2.7) 

Let us recall the relevant properties of the Bargmann transformation. 

1. Ub is a unitary map between L^(R'^) and J- = J-^, the space of all entire holomorphic 
functions f{z) : ^ C such that (here z = x + iy): 

\\f{z)\\l:= f \fiz)\^e-\^\'dxdy = {f,f)F<+oo (2.8) 

where the scalar product {f,g)j^ in J^^ is defined by 

{f,9):F= [ fiz)'^e-\'\"dxdy (2.9) 

Namely, with f{z) := (C/Bn)(z): ||/(z)||^ = \\u{q)h2^^,y 
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2. Let a*, ai be the destruction and creation operators in the variable Xj defined as in 
()1.2p . Let Ni := a*ai be the corresponding number operator, i = 1, . . . ,d. Denote 

d 

N^^^ := ^ A'j the total number operator. Then we have: 



i=l 



UBa*U^^ = z,, UBaiUs^ = —, UbN^Ub^ = Y, Zi— (2.10) 

' i=l * 

so that Hq = N^"^^ + 1. The above operators are defined in their maximal domain in 
Td- Moreover: 

Q{g) := UB{H{g)-l)U^^ = Ub{N^^^ +igalai)Us^ = (2.11) 

9 d . d ^ . 

= zi- Vz2-^ V igz2-^ := Qo + igW 

OZi OZ2 OZi 

defined on the maximal domain. Remark that Spec((5o) = {0, 1, . . . , m, . . .}. The 
eigenvalue Am = "m has multiplicity m + 1. 

3. Let ipk{x) be the normalized eigenvectors of the one-dimensional harmonic oscillator 
in L2(R). Then: 

{UBi^k){z):=ek{z)= J-— z\ fc = 0,l,... (2.12) 

Let now m = 0, 1, 2, . . .. Define: 

fm,h{^l,Z2) ■■= era-h{z2)eh{zi), h = 0,...,m] 
fCm ■■= Span{/,„,h : h = 0,. . . ,m} = 

= Span{eii {z2)ei2 (zi) : + I2 = m} 

Hence the following properties are immediately checked: 

00 

dim/Cm = m + l; /C^ -L /C;, m / /; ^ ICm = ^2 (2-13) 

m=0 

We then have 
Lemma 2.3 

1. For any m = 0,1, . . .: 

Q{g)fm,h = mfm,h + ighfm,h~i^ h = 0,...,m. (2.14) 

2. Let Hm be the orthogonal projection from T2 onto ICm- Then: 
\llrn,Q{g)\ = 0; equivalently, Km reduces Q{g): Q{g)KLrn C Km; 



3. Let Q{g)m ■= Q{g)\Km = 'n.mQ{g)'n.m = ^mQig) = Q{g)'[^m be ICm- component of 

oo 

Q{g). Then Q{g) = @ Q{g)m; 

771=0 

Proof 

1. Just compute the action of Q{g) on fm,h'- 

d d d 
Qi9)fm,h = {zi- h2:27^ \-i9Z2^—)em-h(z2)eh(zi) 

OZi OZ2 OZi 

= {m - h)em~h{z2)eh{zi) + hem~h{z2)eh{zi) (2.15) 
+ ig^h{m- h + l)era-{h~i){z2)eh~i{zi) 
= 'rnfm,h + ig\/ - h + l)fm,h~i 

2. Since the vectors fm,h ■ h = 0, . . . ,m span ICm, by hnearity the above formula entails 
Q[g)K.m c /Cm. 

3. The assertion fohows from 2. above and the completeness relation (|2.13p . 
Proof of Theorem [272] 

We have to prove Assertion 2. 
2. Making /i = in (f2l^ we get: 

Q{9)fm,o = mfmfl, m = 0, 1, . . . 

Hence = m is an eigenvalue of Q{g) with eigenvector fm,o, i.e. with geometric mul- 
tiplicity one. By the unitary equivalence H{g) = U^^{Q{g) + 1)Ub we conclude that 
Am = m + l,m = 0,..., is an eigenvalue of H{g) of geometric multiplicity one, with 
eigenvector U^^fm,o = ''Pm{xi)ipo{x2)- From ()2.14p we read off the matrix representation 
(|2.2p , (j2.3p and we get the statement about the algebraic multiplicity. On account of the 
unitary equivalence ICm = U^^Hm this concludes the proof of the theorem. 

3 A class of non self-adjoint VT symmetric operators with 
real discrete spectrum 

Let Hq be a selfadjoint operator in L^(R'^),d > 1, bounded below (without loss of gener- 
ality, positive) with compact resolvent, and let D{Hq) denote its domain. Let V be the 
parity operator in L^(M'^) defined by 

{ViP){x) = tlj{-x) , VV' G L'^iR'^) , Vx G M'^ . (3.1) 

Let us assume that Hq is P-symmetric, i.e. 

VHi^ = HVi) , G D{Ho) (3.2) 
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and also T-symmetric, i.e. 

(H^){x) = , y^j G D{Ho) , Vx e M'^ . (3.3) 

Let < £i < ^2 < • • • be the increasing sequence of the eigenvalues of TCq. Let rrir denote 
the multiplicity oiij- and ipr,s, -5 = 1, ■■■,'mr, denote linearly independent eigenfunctions 
corresponding to £i , which form a basis of the eigenspace 

A4r '■= SpanjV'r.s : s = 1,...,™,^} (3.4) 

corresponding to £r- 
Definition 3.1 

1. An eigenspace Mr is even (odd) if all basis vectors {ipr,s '■ s = are even 
(odd); i.e., if either Vil)r,s = ipr,s,ys = 1, . . . , rrir, or T'i/Jr,s = — V'r.sj Vs = 1, . . . , m^. 

2. An eigenvalue £r is even (odd) if the corresponding eigenspace Air is even (odd). 

Now, let W G L°°(R'^) be an odd real function, i.e. W{x) = -W{-x), Vx G R"'. Let 
V := iW; clearly V is VT- even, i.e. 

V{-x) = V{x), VxGM'^. (3.5) 

Then, Vg G C, the operator H{g) := Hq + gV defined on D{H{g)) = D{Hq) by 

H{g)i^ = Hoi; + gVi;, e D{Ho) (3.6) 

is closed. More precisely H{g) represents an analytic family of type A of closed operators 
in the sense of Kato ([25J, Ch. VIL2) for 5 G C, with compact resolvents. Thus Spec{H (g)) 
is discrete for all g. For 5 G R the operator H(g) is "PT-symmetric, i.e. 

VH{g)i;{x) = H{g)i^{-x) , W G Z)(/7o) • (3.7) 

Moreover: 

Hig)* = H{-g) (3.8) 
We want to prove the following result. 

Theorem 3.2 Let Hq and W enjoy the above listed properties. Assume furthermore: 

(1) 6:= ^inf(4+i-4) >0; 

(2) Each eigenvalue ir ■ r = 1, . . . is either even or odd. 
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Then if \g\ < ,, ,, — each eigenvalue X{g) of H{g) is real, and thus the spectrum of H{g) 

W*' I loo 

is purely real. 
Example 

The d-dimensional harmonic oscillator with equal frequencies 

d 

2 



k=l 



" ,22 



(3; 



has the properties required by Hq. In this case indeed: 

ir = ujiri + ...+rd + d/2) := uj{r + d/2), = 0, 1, 2, . . . ; k = l,...,d 

with multiplicity = (r + l)*^. Here the corresponding eigenspace is: 

Mr ■■= Span{'4>r,s : s = 1, ...,mr} = Span{ipr^{xi) ■ ■■^nixd) : n + . . . + = r} 

where, as above, ipr{x) is an Hermite function. Now if r is odd the sum r = ri + . . . + r^ 
contains an odd number of odd terms; since "tpsix) is an odd function when s is odd, the 
product i/'ri {xi) ■ ■ ■ tpr^ixd) contains an odd number of odd factors and is therefore odd. ir 
is therefore an odd eigenvalue. An analogous argument shows that ir is an even eigenvalue 
when r is even. Moreover, ir+i — ir = ^ and thus condition (1) above is fulfilled. 

Actually, the above example is a particular case of a more general statement, while for 
d = 2 the above application to the perturbation of harmonic oscillators can be considerably 
strenghtened. 



Theorem 3.3 Let 



2„2 

r- -t- ill 

dx 



2 I ^kXk 



(3.10) 



Assume the frequencies to he rational multiples of a fixed frequency uo > namely: 

cvk = —(^, k = l,...,d (3-11) 

Qk 

where {pk,qk) o,re relatively prime natural numbers. Then: 

(i) If Pk and are both odd, k = 1, . . . ,d, the assumptions of Theorem \3.'^ are fulfilled; 

(ii) If d = 2, the condition \3.11]) pt and qk both odd is also necessary for the validity 
of assumption (2) of Theoren while assumption (1) holds independently of the 
parity of pk, qk- 
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We will now prove Theorem 13.21 in two steps (Propositions [531 and and I3.1(J]) . while the 
proof of Theorem 13.31 is postponed to the next Section. In the first step we show that the 
degenerate Rayleigh-Schrodinger perturbation theory near each eigenvalue Ir is real and 
convergent, with a convergence radius independent of r. Thus there exists p > such that 
all the rrir eigenvalues near ir (counted according to their multiplicity) existing for < /9 
are real for all r. The second step is the proof that H{g) admits no other eigenvalue 
for \g\ < p. To formulate the first step, we recall some relevant notions and results of 
perturbation theory. 

Let go £ Che fixed and let fi be an eigenvalue of H{go). Let c > be sufficiently small 
so that 

Tc = {z ■] z - p\= c] 

encloses no other eigenvalue of H{gQ). Then for \g — g^l small Tc is contained in the 
resolvent set of H{g), p{H{g)) := C \ S^ec{H{g)). Moreover Tc C P, where 

P := {z G C : 36(z) > s.t. {z - H{g))-^ := Rg{z) 
exists and is uniformly bounded for \g — go] < b{z)}. 

Then for \g — go\ sufficiently small 

P{g) = {2m)-^ j Rg{z)dz (3.12) 

is the projection corresponding to the part of the spectrum of H{g) enclosed 

\\Rg{z) - Rg^{z)\\ ^ ^ , as g ^ go (3.13) 

whence 

\\Pig) - Pigo)\\ ^ , asg^go (3.14) 

(see e.g.[25], §VIL1). In particular, if m denotes the multiplicity of p, for g close to H{g) 
has exactly m eigenvalues (counting multiplicity) inside Tc, denoted fisig),s = l,...,m, 
which converge to p as g ^ go- If we denote by A4{g) the range of the projection operator 
P{g), then dimA^(g() = m as g ^ go, and H{g)M.{g) C M{g)- Hence the component 
P{g)H[g)P{g) = P[g)H{g) = H[g)P{g) of H{g) in M{g) has rank m and its eigenvalues 
are precisely ps{g)-,s = 1, m. 

Assume from now on (70 = so that the unperturbed operator is the self-adjoint operator 
Hq := H{0). Let i = £r,r = 1,2,..., be a fixed eigenvalue of Ho, m = rrir its multiplicity 
and ■= '^r,s ■ s = 1, ... ,m be an orthonormal basis in Mr '■= A4r (0). Then there is 
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g{r) > such that the vectors Pr{g)ipr,s ■ s = 1, . . . ,m are a basis in the invariant subspace 
M.r{g) for 1^1 < g{r). We denote 4>r,s{g) : s = 1, . . . ,m the orthonormal basis in M.r{g) 
obtained from Pr{g)ipr,s '■ s = l,...,mr through the Gram- Schmidt orthogonahzation 
procedure. Then the eigenvalues fj,s{g) = ^T,s{g)^s = l,...,m,., are the eigenvalues of the 
rrir X rrir matrix Tr{g) given by: 

{Tr{g))hk ■■= {4>r,h{9), H{g)P{g)4>r^k{g)) = 

{4'r,h{9)^Pr{9)H{g)Pr{g)4>r,k{9)), h, k = 1, . . . , TTlr . 

m 

Let (l)r,s{9) = "lj(9)Pr{g)A,j , a'sjig) G C , s,i = 1, . . . , m^. Then 

m 

{Tr{g))hk = Y alj{g)ali{g){'4}r,j,Pr{-g)H{g)Pr{g)^r,i) , h,k = l,...,m. (3.15) 

Consider now the rrir x ""T-r matrix Br{g) = (-BJi(5))j,z=i,...,m) where 

B';i{g) = {A,j,Pr{-g)H{g)Pr{g)A,i), j, / = 1, . . . , . (3.16) 
Its self-adjointness entails the self-adjointness of Tr{g). We have indeed: 
Lemma 3.4 Let B^ji{g) = BJ^yj, Z = 1, . . . , m,.. Then: 

{Tr{g))hk = {Tr{g))kh, h,k = l,...,mr. 

Proof 

Since Bji{g) = Bfj{g),yj,l we can write: 

rrir rrir 

{TrX9))kh = E '^U9)^UF)B;A9) = E ^I>)"L(5)^Ip(9) 

p,s=l P,s=l 
rrir 

= E al^{g)^j3)B]i{g) = {Tr{g))hk ■ (3.17) 
j,i=i 

and this proves the assertion. 

In other words the selfadjointness of Tr{g)), and thus the reality of the eigenvalues 
ir,s{g) for \g\ < di^), follows from the selfadjointness of Br{g) which will be proved by the 
construction of the Rayleigh-Schodinger perturbation expansion (RSPE) for the operator 
Pr{—g)H{g)Pr{g), which we now briefly recall, following ([25], §11.2.7; here T^^^ = V = 
iW, T^"^ =0„ u>2, D = 0). 

(1) The geometric expansion in powers of g of the resolvent 

oo 

Rgiz) = {z- H{g))-^ = {z-Ho- gV)'^ = ^z) Y^{-gT[V R,{z)T 

n=0 
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is norm convergent for \g\ suitably small. Insertion in (|3.12|) yields the expansion for 

Pig)-- 

oo 

Prig) = E^'^r^"^ P(0) = P.(0) := (3.18) 

n=0 

^i"^ = ^4^/ Ro{z)[VRo{z)rdz, n>l (3.19) 



whence 



Pr{-g)H{g)Pr{g) = £ g^fi:^\ = HoPr (3.20) 

n=0 



where 



fj") = Y,{-lY[PlP^ H^P^'^-P^ + P,(p-iVp,("-p)] , n > 1, P,(-i) = . (3.21) 
and 

= (-1)"+! E V^^'^^)^ _ _ ^ F5^'=")l/5^'="+i) . (3.22) 

fcl + ... + fc„ + l=n 

Here 

5(0) = -P,; 5, = -EP,-/(£,--4); Si''^ = {Sr)\ VA: = 1,2,..., (3.23) 

where Pj is the projection corresponding to the eigenvalue ij of Hq. 

(2) The series (|3.18l3.20p are norm convergent for \g\ < ^^^^^^ — , where dr is the distance 

i = ij. from the rest of the spectrum of Hq. Hence under the present assumptions 
the convergence takes place a fortiori for 

"'-'WU' 

(3) The projection operator Pr{g) is holomorphic for \g\ < p. This entails that its dimen- 
sion is constant throughout the disk. Therefore H{g) admits exactly rrir eigenvalues 
ir^s (counting multiplicities) inside for \g\ < p. 

(4) Hence, for \g\ < p we can write: 

oo 

Brig) = Y.g"gi^\ (g(")),7:=(V'r-,„r^Vr,i), j,l = l,...,mr. (3.25) 

n=0 

We can now formulate the first step: 
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Proposition 3.5 Let ir,r = 1,2,... be an eigenvalue of Hq. Then the nrij. eigenvalues 
(counting multiplicity) Ir^s of H{g) existing for \g\ < p, and converging to ir as g ^ 0, 
are real for \g\ < g{r), 51 E R. 

Proof 

We drop the index r because the argument is r— independent, i.e. we consider the expan- 
sion near the unperturbed eigenvalue i := Ir- Accordingly, we denote hy ips '■= ipr,s the 
corresponding eigenvectors. Let us first consider the case of i even. It is enough to prove 
that = if n is odd and that is selfadjoint (in fact, real symmetric) when n is even. 
These assertions will be proved in Lemma 13.71 and 1 3.91 respectively, which in turn require 
an auxiliary statement. 

Definition 3.6 The product 

U{ki, kn+l) := ^('^i V5('=2)y _ y5(fen)y^{fcn+i) (3.26) 
containing precisely n factors V and n + 1 factors S^^\ j > 0, is called string of length n. 
Then from (j3. 2113. 221) we get: 

n 

= (-1)" T.i-^n(^i%s - igi%s] (3.27) 

where 

iGi%s = {^P„ E U{ki,...,kp+i)Ho J2 n(/ii,...,^„„p+i)^,) (3.28) 

fe;>0 =ii-p;h;>0 

iGS),s = {^„ E n(fci,...,A:p)y E mi,---,K-p+i)^s) (3.29) 

ki + ... + kp 'ii + ... + h,i_p_|.i 

=p-l;fe;>0 =n-p;h;>0 

Now S^''^ is selfadjoint for all k, and V = iW with W{x) G M. Therefore: 

{0i%s = {-ir{ E u{kp+u---,h)i^q,Hon{hi,...,hn-p+i)iJs) 

hi-\'...-\-h^_pj^l—n — p- hi>0 

iQ2!p)<ls = i-ir~'{ E Uikp,...,ki)i;g,VU{hi,...,hn-p+M. 

ki + ... + kp=p-l; fc;>0 
hi + ... + h^_p^l—n — p; hi>0 

Since 5^^^^ ± P, k > 1, in both scalar products (|3.28p and (I3.29P all terms with /ci / or 
hn~p+i / vanish. Hence: 

iGi%s = {-ir{ E U{kp,...,k,)V^Pg,HoU{hi,...,hn^p)V^Ps) (3.30) 

ki + ... + kp—p 
hi+...+hj^_p — n — p 

ki>0,hi>0 

iGi%s = {-ir-\ E U{kp_u---,h)ViPg,VUihi,...,hn-p)V^I;s) (3.31) 

fcl+...+fcp_l— p— 1 

/i]^+...+/i^„p — n — p 
ki>OMi>0 
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We now have: 



Lemma 3.7 Let n be odd, and < p < n. Then, \/ ki, . . . ,kp > 0, V /ii, . . . , hn~p > 0, 
y q, s = 1, . . . ,m: 

{Uikp,...,ki)Vi^g,Hon{hi,...,hn-p)Vi;s) = (3.32) 

(n( Vi, • • • , kiW^Pg, VU{hi, hn-p)Vi^s) = (3.33) 

Proof 

Let us write explicitly (|3.32l I3.33P : 

^^(fcp)y5(fcp-i)y _ _ _ y^^'^i VV'g, HoS^^^W S^^-'W . . . VS^^"-^^Vi;s) = (3.34) 

(5(fcp-i)y^fe-2)y _ ^ ^ vs^'^'^vijq, vs^^^^vs^'^^^v . . . vs^'^^-^^vi^s) = o (3.35) 

Let us now further simplify the notation as follows. We set: 

E ^- = - E (3-36) 

Both series are convergent because \{£j — £)\ > 6 and ^ Pj is convergent. Hence S = 
S-^- © S- and for A; / we have: 

S' = Si^S'_ = {-l)^ E + E ■ (3-37) 

jj^r;£jeven ^ ^ ' jj^r\£jodd ^ ^ ' 

Finally we set 5^"* := S\ := —P. Now, the multiplication by V changes the parity of a 
function, and il'jii^i ^"^^ even. This entails that in both scalar products above S^^'^^ can 
be replaced by S^^^\ S^^^) and so on. The general rule is: S^^^^ can be replaced 

by ( by S^) if and only if j is odd (j is even, respectively). Similarly for the S^^^\ 
Consider first the scalar product in ()3.34p . According to the general rule S'^f*''* coincides 
with S^*"^ if p is even and with S^^'^^ if p is odd. Similarly for S^"^^\ If n is odd p and 
n — p have opposite parity and since Hq does not change the parity of a function the scalar 
product is zero. A similar argument shows that also the scalar product (j3.35p is zero if n 
is odd. Indeed the function in the left hand side has the same parity of the number p—1, 
whereas the function of the right hand side has the same parity of n — p + 1 , and if n is 
odd p — 1 and n — p +1 have opposite parity. This proves the assertion. 

Lemma 3.8 Let n be odd. Then = 0. 

Proof 

It is an immediate consequence of Lemma 13.71 on account of ()3.27|3.30|3.3ip . 
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Lemma 3.9 Let n he even. Then = {Q^'^^)qs for all q, s = 1, . . . ,m. 

Proof 

Once more by ()3.27|3.30|3.3ip we can write for all n (replacing of course V by iW in the 
definition ()3.26p . and denoting IT' the resulting string) 

n 

«"E[ E {-lf{W{ki,...,k,)W^,,H^I\'{hi,...,K-p)W^s) 

p=0 fci + ... + fcp-p;fcj>0 

/i]^+... + 'i7T,_p— p; hj>0 

fci + ... + fcp_]^— p— l;fcj>0 
h-^^--\-...-\-h^_p — n—p; hj>0 

n 

w'Et E {-ir-p{Hou'{h,...,kn-p)w^iJg,u'{hi,...,hp)w^Ps) 

p=0 fci + ... + fc„_p=n-p; 

hi + ... + hp=p:hj>0,kj>0 

^ (-l)"-P+i(W^n'(A:i,...,A;„_p)VFV„n'(/ii),...,Vi)W^V'5)] 

fc]^ + .. . + fc^_p — n— p;fcj- >0 
/il + ... + /ip_l-p-l;/t^>0 

= (3.38) 

To obtain the second equality in (|3.38|) we have used the selfadjointness of Hq and W and 
we have renamed the indices, exchanging p and n — p in the first scalar product, and p—1 
and n — p in the second scalar product. Finally, to obtain the last equality in (I3.38P notice 
that (—1)^ = (—1)""^ since n is even. 
Remarks 

1. It is worth noticing that if the ips, s = 1, . . . ,m, are chosen to be real valued then 
iG^"'^)qs £ RjVj, because W is also real valued and the operators S^''^ map real 
valued functions into real valued functions. 

2. The argument yielding the real nature of the perturbation expansion is independent 
of its convergence, namely it holds for all odd potentials V for which the perturbation 
expansion exists to all orders. In particular, it holds when V is any odd polynomial, 
i.e. for any odd anharmonic oscillators in any dimension d. 

We now proceed to prove that the eigenvalues ir,s{g) are real \/ g £ R, \g\ < p. 

Proposition 3.10 The eigenvalues ir,s, r = 1,2, . . ., s = 1, . . . ,mr are holomorphic for 
\g\ < p and real for g GR, \g\ < p. 

Proof 

The vectors Ur{g)Prfpr,k = Ur{g)tpr,k : k = 1, . . . ,mr represent a basis of M.r{g) for all 
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\g\ < p ([25j. §11.4.2). Here the similarity operator Ur{g)Pr is recursively defined in the 
following way: 

oo 

Ur{g)Pr = Pr + Y. ^r^^ 9^ kU^^^ = kP^,^^ + [k - l)P^^-^^UP + ... P^U^^'^^ (3.39) 

k=l 

We denote Xr,s{9) : s = 1, . . . , rrir the orthonormal basis in M.r{g) obtained from Ur{g)Tpr,s, 
s = 1, . . . ,mr through the Gram-Schmidt orthogonalization procedure. Then the eigen- 
values (.r,s{g)^ s = 1, rrir, are the eigenvalues of the m x m matrix Xr{g) given by: 

{Xr{g))hk ■■= {Xr,h{g),H{g)P{g)xr,k{g)) 

= {Xr,h{g),H{g)xr,k{g)), h,k = l,...,mr. (3.40) 

because P{g)Xr,h{g) = Xr,h-, h = 1, . . . , nir- For \g\ < g{r) the orthonormal vectors Xr,h{g) '■ 
h = 1, . . . , rrir are linear combinations of the orthonormal vectors (pr,hig) ■ h = 1, . . . , rrir 
defined above. Since Xr{g) and Tr{g) represent the same operator on two different or- 
thonormal basis, if either one is self-adjoint the second must enjoy the same property. 
Hence the matrix {Xr{g))hk is self-adjoint, \g\ < g{r), g G R. Expand now {Xr{g))hk in 
power series: 

oo 

{Xr{g))hk = X! {^r,m)hkg^ 

The series converges for \g\ < p. It follows indeed by the standard Gram-Schmidt pro- 
cedure (we omit the details) that it can be written as the quotient of two functions of 
g involving only linear combinations of scalar products of the operators Pr{g) on vec- 
tors independent of g] the denominator never vanishes for \g\ < p by construction, on 
account of the linear independence of the vectors Ur{g)ipr,s, s = 1, . . . ,mr when \g\ < p. 
Now it necessarily follows from the self-adjointness of {Xr{g))hk, valid for \g\ < g[r) that 
{dr,m)hk = i^r,m)kh, m = 0, 1, . . .. Hence the matrix Xr{g) is self-adjoint for \g\ < p, g G R, 
and thus the eigenvalues ir,s ai'e real in the same domain. This proves the assertion. 

Proof of Theorem [Q] 

We have seen that the RSPE associated with the £r-group of eigenvalues ir,sig),s = 
1, . . . ,mr, of H{g) which converge to ir as g ^ 0, have radius of convergence no smaller 
than p. Hence, \/g £ R such that \g\ < p, H[g) admits a sequence of real eigenvalues 
^r,s{g), s = 1, . . . , rrir, r E N. We want to prove that for \g\ < p, g £ R, H{g) has no other 
eigenvalues. Thus all its eigenvalues are real. To this end, for any r G N let Qr denote 
the square centered at Ir with side 25. Then if G R, \g\ < p, and £{g) is an eigenvalue of 
H{g): 

iig) G U Qr. 
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In fact, for any z ^ U^eN Qr we have 

\\gVM^)\\ < l5ll|W^I|oo||i?o(^)|| < p||W^||oo[dist(z,a(i/o))]-' < = 1 (3.41) 

where Roiz) := (Hq - z)~'^ . Thus, z G p{H{g)) and 

R{g, z) ■= {H{g) - z)-^ = Ro{z)[l + gVRo{z)]-^ . 

Now let (7o G be fixed with \g\ < p. Without loss of generality we assume that go > 0. 
Let i{go) be a given eigenvalue of H{go). Then i{go) must be contained in the interior (and 
not on the boundary) of Qn^ for some no G N. Moreover if mo is the multiplicity of ^(50)) 
for g close to go there are mo eigenvalues i^"'^ (5), a = 1. ... , mo, of H{g) which converge to 
i{go) g ^ go and each function £^°'^ {g) represents a branch of one or several holomorphic 
functions which have at most algebraic singularities dX g = go (see [Kato, Thm. VII. 1.8]). 
Let us now follow one of such branches i^°'\g) for < g < go, suppressing the index a 
from now on. First of all we notice that, by continuity, £{g) cannot go out of Qn^ for g 
close to go- Moreover, if we denote the boundary of the square centered at ing with 
side 2t, for < t < 1, we have, for z G T2t and < ^ < ^O; 

\\gVRo{z) II < g[dist{z, a{Ho))]-' < g/t . (3.42) 

Then t > g implies z ^ a{H(g)), i.e. if z G (j{H{g)) H T2t then t < g < go < 1. Hence 
wc observe that as g ^ g^ , i{g) is contained in the square centered at ino and side 2g. 
Suppose that the holomorphic function £(g) is defined on the interval Jgi, 50] with gi > 0. 
We will show that it can be continued up to (7 = 0, and in fact up to = —1. From 
what has been established so far the function f.{g) is bounded as g — >■ gf. Thus, by 
the well known properties on the stability of the eigenvalues of the analytic families of 
operators, £{g) must converge to an eigenvalue £{gi) of H{gi) as g ^ gf and i{gi) is 
contained in the square centered at ino ^^^e 2gi. Repeating the argument starting 
now from i{gi), we can continue £{g) to a holomorphic function on an interval ]52,5i], 
which has at most an algebraic singularity aX g = g2- We build in this way a sequence 
gi > g2 > ■ ■ ■ > gn > ■ ■ ■ which can accumulate only at 5 = — 1. In particular the function 
i{g) is piecewise holomorphic on ] — 1, 1]. But while passing through 5 = 0, £{g) coincides 
with one of the eigenvalues £r,s{g),s = 1, . . . , rrir, generated by an unperturbed eigenvalue 
£r of Hq (namely £no)j which represent rrir real analytic functions defined for 5 g] — 1, 1]. 
Thus, £{go) arises from one of these functions and is therefore real. This concludes the 
proof of the Theorem. 
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4 Perturbation of resonant harmonic oscillators 



Consider again the d-dimensional harmonic oscillator 

1 



dx 



2 2 



2 "^k-^k 



(4.1) 



where now the frequencies uJk > : k = 1, . . . ,d may be different. Theorem 13.31 will be a 
consequence of the following 

Proposition 4.1 The operator j[ j fulfills Assumption (2) of Theorem \ 3.^ if and only 
if the following condition on the frequencies holds: 

(A) 'i k ^ \ {0} such that the components ki : i = 1, . . . ,d have no common divisor, 
and iOiki + . . . + oo^kd = 0, the number 0{k) of ki odd is even. 

Proof 

We first prove the sufficiency part. Let therefore (A) be fulfilled. First recall the obvious 
fact that the rational dependence of the frequency entails the degeneracy of any eigenvalue 
of ()4.ip . In order to show that each eigenvalue 

^ni,..,nd = ^ini + ... UJdUd + ^(Wl + ...UJd) 

of Hq has a definite parity, consider a corresponding eigenfunction 

d 

^E'ni,...,nd(a;i, . . . = n ^ns{Xs) 

s=l 

Now ipns{x) is even or odd according to the parity of Ug, and therefore ^ will be even if 
and only if the number of odd ris is even. Since £ is degenerate, there exist {li, . . . ,ld) / 
(ni, . . . , Ud) such that 

wini + . . . + Udnd = uJih + . . . + Udld =^ (w, k) =0, k := (ni - /i, . . . , - Id) 

and hence the eigenfunction 

d 

'^h,...,la{^l,-..,Xd) = n ^hi^s) 
s=l 

corresponds to the same eigenvalue. The eigenfunctions '^ni,...,na ^-iid ^u, have one 
and the same parity if and only if the number of the odd differences ki is even: in fact, 
an even difference ki = Ui — U does not change the relative parity, while an odd difference 
does. Let us show that if Assumption (A) holds the number of odd differences is even. 
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The case in which ki : i = 1, . . . ,d have no common divisor is the Assumption itself. Let 
therefore ki : i = 1,. . . ,d have a common divisor. If a common divisor is 2, ki is even 
for any i. Hence there are no odd differences. If 2 is not a common divisor, there will 
be an odd common divisor, denoted 6, such that ki = bk'i, where the numbers k'i have no 
common divisor. Now {k',uj) = {k,uj)/b = 0. Hence by the assumptions 0{k') is even. 
Since the multiplication by the odd number b does not change the parity of the k'i, the 
same conclusion applies also to the numbers ki. Thus the total number of odd differences 
does not change after multiplication by b: 0{k) = 0{k') is even. 

Conversely, let us assume that Assumption (A) is violated. Therefore there exists k G 
Z'^ \ {0} such that the numbers ki have no common divisor, {k,u>) = and 0{k) is 
odd. Consider again the eigenfunctions '^ni,...,na{^ij ■ ■ ■ i^rf) ^^'^ • • • ,Xd) cor- 

responding to the same eigenvalue i, with ki = rii — li as above. By construction, the two 
eigenfunctions have opposite parity, and this concludes the proof of the Proposition. 
Proof of Theorem [331 

Let us first prove that Assumption (1) of Theorem 13.21 is fulfilled. Let ii = ^z^,...,;^ and 
(■n = ^ni,...,na denote different eigenvalues. Then, by assumption: 

\en - ill = - h)^ + . . . + (nrf - = 

Qi Qd 

-\{ni - h)piq2 ■■■qd + ■■■ + {nd- ld)Pdqi ■ ■ ■ Qd-ll 



qi---qd 



qi---qd 

Since this lower bound does not depend on the multi-indices (n, /) the assertion is proved. 
Let us now check Assertion (i), namely that if the frequencies have the form tuk = ^Pk/qk 
with and pk odd then Assertion (2) of Theorem (j3.2p holds; namely, all eigenvalues of 
()4.ip have a definite parity. By Proposition 4.2, it is enough to prove that Assumption 
(A) is satisfied. Let indeed (fei, . . . , kd) G Z'' \ {0} be without common divisor and such 
that {oj,k) = 0. Then: 

—ki + ... + —kd = — ^ — {piq2 ■ ■ ■ qdki + ■ ■ ■ + Pdqi ■ ■ ■ qd-ikd) = 
qi qd 91 • • • qd 

■■= — - — (Diki + . . . + Ddkd) = 
qi---qd 

Now the integers : k = 1, . . . ,d are odd; hence the above sum must have an even number 
of terms. The odd terms are those, and only those, containing an odd kf, therefore the 
number of odd ki must be even. Then the result follows by the above Proposition. 
Consider now Assertion (ii) of Theorem 13.31 The only thing left to prove is that the 
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validity of Assumption (A) entails that — = — where di and d2 are odd. Suppose 

indeed — = 7^ where ki is odd and ^2 even, or viceversa. Then ojiki — 002^2 = 0- 
W2 ki 

However this contradicts Assumption (A) which states that the number 0{k) of odd k^ 
must be even. This conchides the proof of Theorem 13.31 

Corollary 4.2 Under the conditions of Theorem \3.3\ on Hq, assume furthermore that the 
matrix {tpr, Wips) ■ r, s = 1, . . . , ttiq is not identically zero for at least one eigenvalue Iq of 
Hq of multiplicity tuq > 1. Then for \g\ < — , H{g) has real eigenvalues if and only 

II ^ lloo 

if p and q are both odd. 
Proof 

The sufficiency part is a particular case of Theorem 13.31 As for the necessity, under the 
present conditions the eigenfunctions have opposite parity. Thefore we can directly apply 
the argument of [TBJ and conclude that if p is even and q odd or viceversa H{g) has a pair 
of complex conjugate eigenvalues near Iq for 5 € K suitably small. 
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